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CHARACTERISTICS  OF  DISCRETE  PROPAGATION  MODES 


ON  A SYSTEM  OF  HORIZONTAL  WIRES 
OVER  A DISSIPATIVE  EARTH 

Steven  W.  Plate,  David  C.  Chang  and  Edward  F.  Kuester 
I . Introduction 

Wave  propagation  along  a system  of  horizontal  wires  parallel  to  the 
earth's  surface  is  of  interest  in  many  areas:  the  transmission  properties 

of  overhead  power  lines,  the  performance  of  antenna  arrays  above  ground, 
the  use  of  Goubau-lines  in  continuous  ground  wave  detection  systems, and 
so  on.  At  present,  most  investigations  of  this  problem  have  been  limited 
to  a low-frequency  approximation  similar  to  that  first  used  for  the  single 
wire  case  by  Carson  (see  the  bibliography  in  [1]).  These  analyses  and 
various  refinements  to  them  [2-14]  are  not  sufficient  to  describe  higher 
frequency  characteristics  of  these  wires.  Only  very  recently  have  more 
rigorous  analyses  of  the  problem  appeared:  in  1959,  Kostenko  [15]  analyzed 
the  case  of  three  identical  conductors  at  the  same  height;  in  1963,  Perel'man 
[16]  analyzed  the  general  m-wire  problem.  A specialization  of  the  latter 
analysis  to  the  two-wire  line  was  done  [17],  as  well  as  an  independent  solu- 
tion of  the  m-wire  problem  [18],  both  in  1965.  Most  recently,  Olsen  and 
Chang  [19]  have  investigated  the  types  of  modes  which  can  exist  on  a two- 
wire  line  over  a conducting  earth  in  light  of  the  earlier  discovery  of  the 
existence  of  so-called  "earth-attached"  modes  on  a single-wire  system  [1,20]. 

In  this  report,  the  equation  for  determining  the  propagation  constants 
of  the  discrete  modes  of  propagation  on  a system  of  m parallel  lines  above 
a dissipative  earth  is  derived.  Our  basic  formulation  of  the  problem  follows 
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closely  that  of  Perel'man  [16].  More  specifically,  we  assume  under  the 
thin-wire  approximation,  that  the  current  on  each  conductor  is  uniformly 
distributed  around  the  surface  of  the  conductor  and  totally  in  the  axial 
direction.  A modal  equation  in  the  form  of  a determinant  is  then  derived 
by  imposing  an  impedance  boundary  condition  in  which  we  assume  that  the 
average  axial  electric  field  equals  the  axial  surface  impedance  multi- 
plied by  the  total  axial  current  on  the  conductor.  Obviously,  this  approxi- 
mation is  good  only  if  the  conductor  radius  is  small  compared  to  the  freespace 
wavelength,  the  distances  between  the  individual  conductors,  and  the  heights 
of  the  conductors  above  the  earth. 

We  find  approximations  to  the  Sommerfeld  integrals  involved  in  the 
modal  equation  that  are  accurate  for  relatively  large  values  of  the  earth's 
refractive  index.  We  then  use  this  appioximate  modal  equation  to  calcu- 
late the  propagation  constants  for  a dual  line  of  bare  wires.  We  find  that 
there  are  two  monofilar  modes  and  one  or  two  bifilar  modes  depending  on  the 
spacing  between  the  wires.  We  also  find  that  a modal  degeneracy  exists  be- 
tween the  two  bifilar  modes  for  a particular  set  of  parameters. 

II.  The  Modal  Equation 

Consider  a system  of  m infinitely  long  parallel  thin  wires  located 
over  a plane  interface  between  two  half  spaces  of  electrical  constants 
Cj,  Uj  and  e2,  y,,  a 2 respectively.  The  j-th  wire  has  a radius  of  a. 
(meters),  and  is  located  at  a height  x=h^  (meters),  and  a position  y=d. 
(meters).  Figure  1 illustrates  the  geometry  of  the  problem.  The  current 
on  the  j-th  wire  is  assumed  to  be  of  the  form  r expfik^z-iut)  where 

a is  the  yet  undetermined  propagation  constant  relative  to  that  of  region 

u 

1,  0)  is  the  angular  frequency  of  the  fields,  and  ^ l el ) 


i 
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Fol lowing  Wait  [21],  expressions  for  the  electric  and  magnetic  fields 
in  the  region  x>0  due  to  a filament  of  current  I located  at  x=x'  , y=y' 
can  be  written  as 


Eu  • kf‘2"l 


• "l:  ■ kf«2"i 


3n, 


an. 


Eix  = ikia  ax-  + 


31^* 


3fli 


; II.  = ik.a- - ie_w- — 

lx  1 3x  1 3y 


an. 


*1  3IT1 

E,  - ik.a- ip,a)- — 

ly  1 3y  1 3x 


311, 


; H,  = ik.a-; 
ly  1 3y 


3n  j 

+ ie  .to- — 
1 3x 


(1) 


The  scalar  potentials  and  n^*  are  the  z components  of  the  electric 

and  magnetic  Hertz  vectors,  respectively.  By  assuming  that  = p,  = p() 
expressions  for  these  potentials  are  found  to  be 


-n.I  , . , 

= -fiq  exp(ikjaz)  {H0UJ  (k^  ((x-x')2  + (y-y')2)  "2  1 


H0(1)  [kjt  ((x  + x’j2  ♦ (y  - y')2  ) * ] 


li^-J  [ 5 ]exP  ['uiki(x  * *’>  - - y')i  dx } (2) 

^ -oo  1 Z 1 


V = l^r^})  exp(ikl0cZ)  | 


T 


Xexpl-Ujkj (x+x’ ) - iXkjCy-y')] 
(u2  + n^UjJCUj  + u2) 


dX  (3) 


Ul  = (X2  - c2)4 
U2  = ‘ 


; -7t/2  < arg  Uj  < tt/2 
; -it/2  < arg  u2  < tt/2 


where 
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C = (1  - ot^)  ; 0 < arg  c < tt 

C = (n2  - a2)^  ; 0 S arg  S < tt 

n n 

The  arguments  of  these  variables  are  defined  so  that  the  electromagnetic 

fields  are  bounded  at  infinity.  In  these  expressions  represents 

the  Ilankel  function  of  the  first  kind  and  of  order  zero,  and  n represents 

the  relative  refractive  index  of  region  2 compared  to  region  1,  defined  as 

n2  = k2  /k2  = e_,(l  ♦ i6)/Cj  with  5 (loss  tangent)  given  as  o^/fb^) 

u 

and  0 <,  arg  n < tt/4  . Finally,  = (u^/e  )^  is  the  intrinsic  impedance 
of  region  1.  We  note  that  the  first  and  second  terms  in  (2)  are  the  direct 
contributions  due  to  the  current  filament  and  its  image  in  region  2,  respec- 
tively. The  integral  in  (2)  and  the  magnetic  potential  given  in  (3)  are  due 
to  the  finite  conductivity  of  region  2. 

With  the  use  of  the  addition  theorem  of  cylindrical  functions,  we  there- 
fore obtain  the  expression  for  the  electric  potential  due  to  a uniform  current 
distributed  on  the  surface  of  the  j -th  wire  as 

"lj  = 4^  exp(iklaZ)Jo(i;ajkl)  fc2Ho1)lCkl((x'hj)2+  (X 

-C2Hq(1)  ( ckj  ( (x  ♦ h.)2  ♦ (y  - d.)2  ) ’s  ] + P [ ajk^x  ♦ h.J.k^y  - d. )] 
-QKkjtx  ♦ hj ) .kj  (y  - d.)]}  (4) 

where 


P(a;X,Y)  = T 


ITT 


C°°  exp(-XUj  - iYA) 
lll  + U2 


dA 


(5) 


2 2 f“  exp (-Xu.  - iYA) 
Q(a;X,Y)  = ^ dA 

’ ’ ’ in  J u2  + n?u 


(6) 
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III.  PROPAGATION  MODES  ON  A DUAL- LINE 

In  the  special  'ase  of  two  identical  bare  wires  of  equal  height  (H ^ , 

the  modal  equation  can  be  factored  to  obtain  two  solutions:  the  monofilar 

one  where  1^  = 1^,  and  the  bifilar  one  where  1^  = -I^.  The  two  independent 
modal  equations  are 

0 = M+(a)  = Jg(CA){^2[H(51)(a)/J0(a)  ' Hq1  3 (?2H)  j ♦ P(a;2H,0) -Q(a;2H,0)} 

± (CA){?2tH^1)(?D)  -hJ1}OT]  + P(a;2H,D)-Q(a;2HfD)}  (12) 

where 

1 

A = Aj  = A2,  H = Hj  = H2,  D = | Dx  - D2 1 , R = f(2H)2  ♦ D2  1 

We  note  that  the  Bessel  function  Jq(CA)  is  only  important  in  the  behavior 
of  M(a)  for  large  £ . For  the  purpose  of  finding  the  roots  of  M(a), 
the  term  J (?A)  may  be  approximated  by  unity  under  the  thin  wire  assumption 
that  A « 1. 

The  method  we  adopt  in  this  report  to  find  the  zeros  of  M+ (a)  is  an 
iterative  scheme  similar  to  the  well-known  Newton's  method.  To  use  this  method 
the  derivative  of  M(a)  with  respect  to  a is  needed.  This  is  given  by 

M±'(a)  = {-2a[H^1)(i;A)-H^1)(^2H)]  * a? [AH^1 3 (r,A)  - 2H}1 } (C2H)  ] 

+ P'(a;2H,0)  - Q'(a;2H,0) 

± {-2afH^1)(CD)  -H^^CCR)]  + a^DH^1 3 (CD)  - RhJ1  5 [£R) ] 

+ P^at; 2H,D)  - Q'(a;2H,D))  (13) 


where 
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P1 (a;X,Y) 


-2a  r 1 + XU2 


t — I j- r-  exp(-u.X  - iAY)dA 

ITT  J U U (U,+U,)  1 J 

- oo  1 L 1 Z 


(14) 


Q'(a;X,Y)  = 


2a  7 r.  J>„  °t2(ul+n2u2)  ,«xp(-u  X-iAY 

j [2u^  - a X , | ^ dA 


ITT 


2 J 2 
u2 (u2+n  u2)  Uj (u2  + n Uj) 


(15) 


The  (j+l)-th  approximation  to  the  root  a^+J  is  then  obtained  from  the 
previous  approximation  a^  using  the  equation 

“j  + l = “j  ' M± C<X j J /M± (0t j 5 * “ 0)  (16) 

where  a(1  is  an  initial  guess  to  the  root,  and  the  iteration  is  stopped  when 
|M(a^)|  < e (where  e is  the  desired  accuracy  of  the  root).  Obviously  the 
iteration  converges  only  if  the  initial  guess  is  close  enough  to  the  actual 
root.  The  method  used  in  this  report  to  locate  the  roots  is  to  find  one  for 
a particular  set  of  parameters  (by  trial  and  error  or  from  previously  known  results 
and  then  vary  the  parameters  by  small  increments,  finding  the  root  at  each  step 
until  the  desired  set  of  parameters  is  reached.  At  each  step  the  initial  guess 
a^  is  the  root  of  the  previous  set  of  parameters. 

The  integrals  P and  Q and  their  derivatives  as  given  in  (5),  (6),  (14), 
and  (15),  in  principle  may  be  evaluated  numerically,  however  it  is  more 
efficient  to  find  approximate  expressions  initially  in  order  to  avoid  excessive 
computation  in  the  complex  a-plane.  In  the  following  section,  we  shall  discuss 
various  approximations  valid  for  different  ranges  of  parameters  involved. 
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IV.  APPROXIMATIONS  FOR  THE  INTEGRALS 
P(a;X,Y)  and  Q(a;X,Y) 

Using  methods  similar  to  those  employed  by  Olsen  and  Chang  [20],  the 
integrals  P and  Q can  be  approximated  in  terms  of  known  functions  for 
the  common  case  where  the  wire  height  above  the  ground  is  greater  than  the 
skin-depth  of  the  ground,  that  is,  |n|  II  » 1.  The  integral  P as  given 
in  (5)  can  be  written  as 

OO 

2 r 


P(a;X,Y)  = - j^2  (u1-u2)exp(-UjX  - iXY)dX 


(17) 


where 


N = (n2  - 1)^  ; 0 <.  arg  N < ir/4 


In  most  cases,  we  expect  the  useful  solutions  to  the  modal  equation 
M±(a)  =0  to  be  located  near  a = 1.  This  means  that  C is  small  so  the 
integrand  of  (17)  decays  as  exp(-X|A|)  away  from  the  point  X = 0.  The 
major  contribution  to  the  integral  is  then  from  the  region  X =»  o.  We 
therefore  expand  u^  in  a Taylor  series  about  X =0  to  obtain 

...  iX2  iX4 
u2  = -lCn  + 2C  + — + 


n 8C 


(18) 


n 


This  series  converges  only  when  |X|  <1^1,  however  the  integrand  of  (17) 
will  have  decayed  by  the  factor  exp(-x|c  |)  outside  these  limits.  Therefore 
an  approximation  based  on  the  first  few  terms  of  (18)  appears  to  be  valid  for 
|Snl  >>  I Sl  an<*  kJX  » 1.  By  using  the  first  term  of  (18),  we  may 
rewrite  (17)  as 


iirN  -a, 


P(a;X,Y)  = —~  f [u  + ir  ]exp(-u.X  - iXY)dA 

/in  i 

OO 

f [-i?n  - u2]exp(- ijX  - iXY)dX 


iTTN 

P„(a;X,Y)  + e 


(19) 


P0 
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Neglecting  £p(),we  obtain  a first  order  approximation  to  P as 


P(a;X,Y)  - P0(u;X,Y) 


= { [32/3X2  - i£  3/3X]  f fexpf-UjX  - iAY)/u  ]dA 

iirtv  n J ™ 


= (2/N2){cH1(1)(?R)tiCnX/R  ♦ (X2-Y2)/R3] 


where  R 


- Crx2/R2)H^1)(CR)  (20) 

2 2 ^ 

(X  + Y ) . In  deriving  (19)  the  following  identity  has  been  used. 


H^fcR)  = 1 ^^[expC-UjX  - iAY)/uJ]dX  (21 

An  upper  bound  for  the  absolute  error  in  this  approximation  is  found  in 
Appendix  A to  be 


lePi)l  < 


4[2*26X  62X2  ♦ 63X5/3] 

Hn\|x3 


where 


|0  (Re  z < 0) 

{ j i 2 

i (Re  Z)  (Re  z > 0) 


From  this  expression  it  is  apparent  that  the  error  is  small  if 
| n3 1 X3  » 8/tt  and  if  |n3|  »463/(3tt)  . For  example  say  n = 5.3  + i.95 
and  Re  z2  < -04,  then  the  error  for  h = .lA  is  less  than  9 x 10~  and 
for  h=  . 5\  the  error  is  less  than  7 x 10  3 . These  error  estimates  are 
conservative  and  as  will  be  shown  in  a later  section,  the  results  obtained 
using  this  approximation  are  generally  much  better  than  what  is  indicated  by 
the  error  bound.  For  the  interested  reader,  we  have  also  included  in 
Appendix  A higher  order  approximations  and  their  error  estimates. 


w 


-11- 

Following  a similar  procedure,  we  now  rearrange  the  integral  Q as 

OO  2 

2 2 f u -u./n 

Q(a;X,Y)  = j -----  %■  exp(-u.  X - iXY)dX 


iir(n  -1)  00  X -X, 


(23) 


where 


2 a?  i 

Xp  = (C  -1/n  ) ; 0 < arg  Xp  < tr 


n = (n*-  ♦ 1)^  ; 0 < 


< arg  n < 7r/4  . 


Since  the  integrand  of  (23)  has  a pair  of  poles  at  ±X^,  it  is  more  convenient 
to  expand  u 2 around  X = X instead  of  around  X = o,  so  that 


2 2 
• r \ *•  \ *■ 


2 2 2 ? 


•in 

u2  = ~TT 


-in2  U*  ‘On  i(x  -X„)  » 


2n 


8n 


(24) 


Based  upon  (24) , we  may  rewrite  (23)  as 


2oc2n2  f ul+i/n  r -i/n-u  /n- 

Q(a;X,Y)  = - ■ 7 { / -5 =-  exp(-u,X  - iXY)dX  ♦ / = exp(-UlX-iXY)dX } 


u.+i/n 


iir(n  -1)  “ X - X 


-i/n  - u2/n 

wr  ~72  71 

00  X - X 


3 Q0(«;X,Y)  ♦ Eqo 


(25) 


where 


Q0(a;X,Y) 


2a2n2 


OC 

I 


exp(-u^X  - iXY) 


iir(n4-l)  J-  V17^ 


dX 


and 


2a 


'Qo 


in (n  -1) 


/ 


exp(-UjX  - iXY) 

7~T7 

u2~in  /n 


dX 


We  note  that  the  integrand  of  Qq  contains  a pair  of  poles  at  X = ± X , 

whereas  the  integrand  of  e_  does  not.  According  to  Olsen  and  Chang  [20] 

y° 

2 a 2 

this  pole  will  cause  a singularity  in  Q where  X^  = 0 (i.e.  at  a = 1 -1/n  ). 
Since  is  small  compared  to  Qq  and  has  a nonsingular  integrand,  we 

may  neglect  it  to  obtain  the  first  order  approximation 


-n- 


Q(a;X,Y)  * Qo(a;X,Y) 


,2  2 

2a  n 


ifl(n  -1)  00 


f ^ /"exp(-u  X-iXY) 

{ J [exp(-UjX  - iXY) /u^JdX  + i/n  J u (u  i/fi)  dX  1 


2slt—  H^fcR)  ♦ W(a;X,Y)  (26a) 

(n  -1)  u ir(n  -l)n 


where  W(a;X,Y)  can  be  written  as 


W(a;X,Y)  = exp( 


00 

■iX/n)  {J 


[exp(-X(Uj-  i/n))-l] 


UjCUj-i/n) 


exp(-iXY)d  X 


♦ / exP(~iXY)  ■ ■ dX  } 

u^Uj-i/n) 


= exp(-iX/n){Wx(a;X,Y)  ♦ WQ(a;Y)} 


(26b) 


We  note  that  the  integrand  of  Wx  has  a pair  of  branch  cuts  in  the 
complex  X-plane  due  to  the  definition  of  Uj  but  does  not  have  any  poles. 
According  to  Olsen  and  Chang  [20]  this  results  in  a pair  of  branch  cuts  in  the 
a-plane  due  to  the  motion  of  the  branch  points  of  Uj  (at  X = ± £)  crossing  the 
real  axis  in  the  X plane.  The  branch  points  in  the  a-plane  are  at  a = ± 1 
and  the  cut  is  defined  by  those  points  where  c is  a real  number. 

The  integrand  of  WQ  has,  in  addition,  two  poles  at  X = tX^.  Again  we 
know  from  the  work  of  Olsen  and  Chang  [20]  that  these  poles  can  cause  branch 
cuts  in  the  a-plane  whenever  they  cross  the  real  axis  in  the  X-plane  because 
of  the  discontinuity  in  the  residue  calculation  at  t A . These  branch  cuts 

2 a2  5 P 

in  the  a-plane  are  located  at  a * t (1  - S - 1/fi  ) where  S is  any  real 


number. 
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An  expression  for  WQ(a;Y)  is  derived  in  Appendix  B.  One  form  of  the 
expression  is 

WQ(a;Y)  = (2/A  )cos  (A  Y)  [fcnO; ) -Hn(l/n  - iA  ) ] - (tt ,'A  )sin(A  Y) 


I 

tt/(A  n)  f sin[A  (Y-s)]H^1^(i;s'ds 
P Jr i P u 


(27) 


where  the  principal  values  of  the  logarithms  are  chosen.  A series  expansion 
for  the  finite  integral  used  above  is  given  as 

“ ^nY)m 

(iY/2)  l b — [exp(-iA  Y)  - (-1)  exp(iA  Y)]  I (c Y ) (28) 

m=0  p pm 


I 

f 

&■ 


I 

I 


where  Im(^Y)  is  expressed  in  terms  of  known  functions  in  equations  (B.14), 
(B. 15) , and  (B.16). 

An  expressianfor  Wx(a;X,Y)  is  found  in  Appendix  C to  be 


Wx(a;X,Y)  = 


* . 

-iir  £ expfis/nJHp1^  [^(s2+Y2)!]ds 


(29) 


A series  representation  of  this  integral  can  be  shown  to  be 


" fm(C.X2/2*Y2)x2m+1 

W (a;X,Y)  = l -2 


m=0 


ml 


^(iX/n) 


(30) 


where  f^  and  I are  again  expressed  in  terms  of  known  functions  in 
Appendix  C.  We  note  that  such  a series  does  not  converge  well  when  Y small 
compared  to  X.  In  order  to  find  an  expression  that  is  good  in  this  region, 
W(a;X,Y)  can  be  rearranged  as 


W(a 


O 

;x,Y)  = f 


exp(-UjX)cos (AY) 


dA 


oo  Uj  (Uj-i/n) 


(31) 


-I 


An  expression  for  (31)  is  found  in  Appendix  D by  expanding  the  cosine 
term  into  a Taylor  series,  and  integrating  this  series  term  by  term.  The 
result  is  repeated  here  as 

00  ( 1lmY2,n 

W (a ; X , Y ) = cos(A  Y)W(ct;X,0)  - iu  \ — R (a;X)  (32) 

p m=l  (2m)!  m 


where 


W(a;X,0) 


exp(-iX/n)  { -iir  I exp (is/n)!!^1  •*  (Cs)ds 
-'0 


♦ (2/y  [*n(0  -fcn(l/n  -iA  )]}. 

Explicit  expressions  for  W(a;X,0)  can  be  found  in  Appendix  B and  will  not  be 


repeated  here.  The  terms  Rm  are  expressed  in  terms  of  known  functions  in 
(D.4),  and  a series  expansion  for  the  finite  integral  is  given  as 


00  v /~>m 

y ]—  * i (a) 

m=0  m! 


where  Im(£X)  is  again  given  in  (B.14),  (B.15),  and  (B.16).  It  should  be 

noted  that  the  series  in  (32)  converges  only  if  Y<X,  so  in  computing  W 
2 2 

we  use  (32)  if  X > 2Y  , otherwise  we  use  (27),  (28),  and  (30)  inserted 


into  (26b). 

In  the  preceding  approximations  on  the  integral  Q,  it  is  assumed  that 
in  summing  each  of  the  series  enough  terms  can  be  included  to  obtain  any 
desired  accuracy.  As  evident  from  (25),  the  term  is  omitted  in  the 

derivation.  An  upper  bound  for  the  error  due  to  neglecting  £_  is  found 
in  Appendix  E to  be 

> . . 4 a2n  (1  +6X) 

I V < n 71“77~2  — x — 
x (n  -l)n 


0 (Re  C < 0)  , 

(Re  cV  (Re  t2  > 0)  . 


where 


-15- 


From  this  expression  it  is  apparent  that  the  error  is  small  if 

r r 

|n‘>|x  » 4/tt  and  |n°|  » 46/tt . For  example  say  n=5.3  + i.95  and 

2 -4 

Real  6 < .04,  then  the  error  for  h=0.lX  is  less  than  3 x 10  and  for 

h=0.5A  the  error  is  less  than  6x10  We  also  note  that  Q is  of  the 

2 3 

, order  1/n  , so  the  relative  error  is  of  the  order  1/n  . This  error 

estimate  is  again  conservative  and  the  results  we  obtain  in  the  following 

section  will  show  that  the  error  is  usually  less  than  indicated  by  these 

estimates.  Similar  to  the  evaluation  of  P(a)  we  have  also  included  in 

Appendix  E higher  order  approximations  and  their  error  estimates.  It  should 

be  noted  the  series  for  P and  Q,  generated  by  including  the  higher  order 

terms  of  u^  are  only  asymptotic,  so  in  general  one  does  not  necessarily 

improve  the  accuracy  by  including  more  terms. 

I i 
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V.  NUMERICAL  RESULTS 

We  have  developed  a computer  program  to  compute  the  roots  of  the  dual- 
line modal  equation  (12)  using  the  first  order  approximations  to  the  Sommerfeld 
integrals  P(a;X,Y)  and  Q(<*;X,Y).  To  test  the  accuracy  of  these  approxi- 
mations we  compared  the  values  of  the  propagation  constants  found  using  the 
first  order  approximations,  to  the  values  found  using  a numerical  integration 

of  the  Sommerfeld  integrals.  Typical  results  given  in  Table  1 show  that  the 

-4 

accuracy  of  the  approximations  is  quite  good  (on  the  order  of  10  or  less), 
even  for  the  case  of  a poorly  conducting  earth  with  |n|  =5.4  considered 
here.  It  is  assumed  that  subsequent  values  obtained  using  the  approximate 
modal  equation  are  of  the  same  order  of  accuracy  as  the  values  in  Table  1. 

In  Fig.  2 the  roots  of  the  modal  equation  are  plotted  for  several 
values  of  the  wire  spacing  d,  with  the  wire  height  h=  0.2X,  radius 
a = . 005A  , and  refractive  index  n = 5.3  + i0.95.  As  expected  for  large 
spacings  there  are  monofilar  and bi filar  modes  with  propagation  constants 
close  to  the  values  of  the  single  wire  modes.  As  the  spacing  is  decreased 
the  attenuation  of  the  "quasi  TEM"  monofilar  mode  increased  until  the  spacing 
is  approximately  equal  to  the  wire  height.  After  this  point  the  propagation 
constant  approaches  that  of  a single  wire  of  radius  equal  to  the  geometric 
mean  of  a and  d.  This  can  be  shown  directly  from  the  modal  equation. 

The  "earth-attached"  monofilar  mode  is  relatively  insensitive  to  the  spacing 
of  the  wires.  The  "quasi  TEM"  bifilar  mode  is  affected  less  by  the  earth 
as  the  spacing  decreases,  because  the  fields  are  concentrated  between  the 
two  wires.  It  should  be  noted  that  the  "earth-attached"  monofilar  mode 
is  less  attenuated  than  the  "quasi  TEM"  bifilar  mode  for  spacings  larger  than 


Mode  trajectories  as  a function  of  wire  height  for 
various  wire  spacings 
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the  wire  height;  however  for  very  small  spacings  the  bi filar  mode  has  the 
lower  attenuation.  The  fourth  mode,  an  "earth-attached"  bifilar  mode,  can 
exist  for  large  wire  spacings  only.  This  mode  disappears  into  the  improper 
Riemann  sheet  at  spacings  of  1.5  to  2.0  wavelengths  for  wire  heights  of  0.1 
to  0.4  wavelengths.  The  disappearance  of  this  mode  is  due  to  the  cancelling 
of  the  singular  portions  of  Q(a;2H,0)  and  Q(a;2H,D)  in  the  modal  equation 
(12).  Note  that  it  is  possible  for  the  "earth-attached"  bifilar  mode  to  have 
a lower  attenuation  than  the  "quasi  TEM"  bifilar  mode. 

Figure  3 shows  a similar  plot  of  the  roots  of  the  modal  equation  for  a 
height  of  h=0.3A.  The  movement  of  the  monofilar  mode  is  similar  to  that 
of  Fig.  2;  however  the  bifilar  mode  which  disappears  and  the  one  which 
becomes  TEM  as  the  spacing  is  decreased  seem  to  have  been  interchanged.  This 
indicates  that  there  must  be  a degeneracy  between  the  two  bifilar  modes  at 
some  height  between  0.15A  and  0.3A  . 

Figure  4 of  the  bifilar  modes  is  a plot  for  several  heights  between 
l). 15A  and  0. 3A  and  for  various  wire  spacings.  This  figure  shows  that  the 
degeneracy  occurs  at  a wire  height  between  0.25A  and  0.3A  and  at  a wire 
spacing  of  about  2.5A  . This  degeneracy  makes  it  difficult  to  label  the 
modes  as  being  either  "earth-attached"  or  "quasi  TEM",  because  these  modes 
can  be  transformed  continuously  into  each  other  by  varying  the  spacing  and 
height  of  the  wires. 

Figure  5 is  a plot  of  the  modes  for  changing  heights  at  several  fixed 
•spacings.  This  shows  that  the  bifilar  mode  that  exists  for  small  spacings 
is  transformed  into  a TEM  type  mode  as  the  height  is  increased.  One  set  of 
monofilar  modes  moves  from  the  branch  point  at  afip  = n/n  to  a = 1 as  the 
height  is  increased.  The  other  set  of  monof i lar  modes  moves  toward  the  branch 
point  as  the  height  is  increased. 
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VI.  CONCLUSION 

In  this  report,  we  have  investigated  the  modes  of  propagation  along  a 
two-wire  line  parallel  to  and  above  the  surface  of  a finitely  conducting 
earth.  Due  to  interaction  between  the  two  wires  as  well  as  between  the 
wires  and  the  earth,  the  mode  structure  is  more  complex  than  that  which 
would  be  found  in  the  case  of  a perfectly  conducting  ground  - a single 
monofilar  and  a single  bifilar  mode.  Moreover,  the  existence  of  modal 
degeneration,  similar  to  that  discovered  in  [1]  for  the  single-wire  line, 
has  been  demonstrated. 

Systematic  analytic  approximations  have  been  derived  for  the  Sommerfeld 
integrals  which  enter  into  the  modal  equation,  as  well  as  rigorous  error 
bounds  for  them.  These  expressions  allow  a great  savings  in  the  computer 
time  required  to  numerically  determine  the  modal  propagation  constants. 


Acknowledgments 


We  wish  to  thank  R.G.  Olsen  of  Washington  State  University,  and 
W.  Rotman  of  RADC/ET  for  many  interesting  discussions  and  suggestions. 


-24- 

REFERENCES 


[1]  E.  F.  Kuester  and  D.  C.  Chang,  "Modal  representation  of  a horizontal 

wire  above  a finitely  conducting  earth,"  Tech.  Rept.  No.  1,  (RADC-TR- 76-287) 
Dept,  of  Elec.  Eng.,  Univ.  of  Colo.,  Boulder,  Colo.  (1976). 

[2]  R.  M.  Foster,  "Mutual  impedance  of  grounded  wires  lying  on  the  surface 
of  the  earth."  Bell  Syst.  Tech.  J.  v.  10,  pp.  408-419  (1931). 

[3]  M.  C.  Gray,  "Mutual  impedance  of  grounded  wires  lying  on  the  surface  of 
the  earth  when  the  conductivity  varies  exponentially  with  depth,"  Physics 
v.  5,  pp.  35-37  (1934). 

[4]  A.  S.  Knyazev,  "A  practical  method  of  calculating  the  resistance  of 
line  conductors  with  the  effect  of  the  real  earth  taken  into  account," 
Radiotekhnika  v.  15,  no.  9,  pp.  21-32  (1960).  [In  Russian;  Eng.  transl. 
in  Radio  Engineering  v.  15,  no.  9,  pp.  25-41  (I960)]. 

[5]  M.  V.  Kostenko,  "Propagation  of  electromagnetic  waves  along  multiconductor 
lines,"  Elektrichestvo,  no.  11,  pp.  8-12  (1960)  [In  Russian]. 

[6]  S.  P.  Belousov  and  V.  G.  Yampol'skii,  "A  twin-wire  traveling-wave  antenna," 
Elektrosvyaz' , no.  5,  pp.  24-30  (1962)  [In  Russian;  Eng.  transl.  in  Tele 
communications,  no.  5,  pp.  23-29  (1962)]  . 

[7]  M.  V.  Kostenko  and  L.  S.  Perel'man,  "On  the  calculation  of  wave  processes 
on  multiconductor  lines,"  Izv.  Akad.  Nauk  SSSR — Energetika  i Transport, 
no.  6,  pp.  721-730  (1963)  [In  Russian]. 

[8]  G.  A.  Lavrov  and  A.  S.  Knyazev,  Prizemnye  i Podzemnye  Antenny,  Moscow: 
Sovetskoe  Radio,  1965  [In  Russian] . 


[9]  L.  M.  Wedepohl  and  R.  G.  Wasley,  "Wave  propagation  in  multiconductor 
overhead  lines,"  Proc.  IEE  (London),  v.  113,  pp.  627-632  (1966). 

[101  M.  Krakowski,  "Active  power  loss  in  the  ferromagnetic  medium  due  to 
currents  in  multiconductor  systems,"  Archlwum  Elektrotechniki,  v.  19, 
pp.  85-105  (1970). 

[11]  M.  Krakowski,  "Distribution  of  active  power  on  surface  of  ferromagnetic 
medium  due  to  currents  in  parallel  conductors,"  Archiwum  Elektrotechniki 
v.  19,  pp.  525-538  (1970). 

[12]  M.  Krakowski,  "Shielding  effect  in  a system  due  to  currents  in  parallel 
conductors,"  Archlwum  Elektrotechniki  v.  21,  pp.  107-123  (1972). 

[13]  M.  Krakowski,  "Eddy  current  losses,  in  a non-ferromagnetic  plate  due  to 
currents  in  a multiconductor  system,"  Archlwum  Elektrotechniki  v.  22, 
pp.  321-350  (1973). 


4 


.v> 


-25- 

[14]  M.  C.  Perz  and  M.  R.  Raghuveer,  "Generalized  derivation  of  fields, 
and  impedance  correction  factors  of  lossy  transmission  lines.  Part  II. 

Lossy  conductors  above  lossy  ground,"  IEEE  Trans.  Power  App.  Syst.  v.  93, 
pp.  1832-1841  (1974). 

[15]  M.  V.  Kostenko,  "Propagation  of  sinusoidal  oscillations  along  a three- 
wire  line  with  horizontal  wire  spacing,"  Elektrichestvo , no.  8,  pp.  1-8 
(1959)  [In  Russian;  also  available  in  Czech  in  Elektrotechn.  Obzor  v.  50, 
pp.  11-18  (1961)] 

[16]  L.  S.  Perel'man,  "Details  of  the  theory  of  wave  propagation  along  multi- 
conductor transmission  lines  in  connection  with  some  engineering  problems," 
Izv.  Nauchno-lssled.  Inst.  Postoyannogo  To ka  (Leningrad),  no.  10,  pp.  103-120 
(1963)  [In  Russian]. 

[17]  F.  G.  Kaidanov,  et  al.,  "Precise  determination  of  transmission-line 
constants  and  analysis  of  the  telegraph  equations  for  the  example  of 
a two-conductor  transmission  line,"  Elektrichestvo,  no.  3,  pp.  15-21 
(1965)  [in  Russian;  Eng.  transl.  in  Elec.  Technol.  USSR,  v.  2,  pn.  179- 
198  (1965)] 

[18]  I.  Vago,  "Over-ground  power  transmission  line  systems,"  Periodica  Poly- 
technics Elec.  EnR.  v.  9,  pp.  67-96  (1965)  [from  orig.  article  in 
Hungarian  in  Elektrotechnika  v.  58,  pp.  190-194  and  236-244  (1965)]. 

[19]  R.  G.  Olsen  and  D.  C.  Chang,  "Propagation  of  modes  on  a pair  of  infinite 
wires  above  a conducting  earth,"  1975  Internatl.  AP-S  Symposium,  pp.  427-430 
(1975). 

[20]  R.  G.  Olsen  and  D.  C.  Chang,  "Electromagnetic  characteristics  of  a hori- 
zontal wire  above  a dissipative  earth — Part  I:  Propagation  of  transmission- 

line and  fast-wave  modes,"  Sci.  Rept.  No.  3,  (NOAA-N22-126-72)  Dept,  of 
Elec.  Eng.,  Univ.  of  Colo.,  Boulder,  Colo.  (1973). 

[ 21 ] J.  R.  Wait,  " Theory  of  wave  propagation  along  a thin  wire  parallel  to 
an  interface,"  Radio  Science  v.  7,  pp.  675-679  (1972). 


-26- 


APPENDIX  A 

In  this  Appendix  the  higher  order  terms  of  P are  derived  and  error 
bounds  for  these  are  found.  An  (M+l)th  order  approximation  is  obtained 

by  extracting  the  first  M+l  terms  in  the  expansion  of  u„  in  (18)  from  u„ 

l i 

in  (17).  P can  then  be  written  as 


M 

P(a;X;Y)  = l p (a;X;Y)  + eDM 
, n k PM 

k=0 


(A.l) 


where  Pq  is  given  in  (20),  and  for  k>0 


2U 

Pk(a;X;Y)  = 


1 sAv**10 


and  where 


and 


N k! 


Ck  = (‘  2 ) • • •<! " k) 


CD 

IR(a;X;Y)  = (iiT)-1  J X2k  exp(-u1X-iXY)dA 


(A. 2) 


(A.  3) 


(iTT)"1{[-32/9Y2]k  [-3/3X]  J [exp(-u1X-iAY)/u1  ] dA } 


= {[-32/3Y2]k  a(X2+Y2)-!s  H(J)fC(X2+Y2)ljJ}  (A. 4) 

with  the  use  of  the  identity  given  in  (21).  In  particular,  for  the 
2- terms  expansion 

P. (a;X;Y)  = (i/N2^  ) { (CX/R5)H(J) <£R)  [6Y2-2X2-£2Y2R2]  + (C2X2/R4)H(1) (<;R) [X2-3Y2] } 
l n l o 

(A.  5) 

Combining  (A. 5)  and  (19)  we  obtain  the  following  expression  for  P(a;X;Y): 

P(a;X;Y)  - \ — { -4  H(}}  (r,R)  [X2-Y2+U  XR2  + 1X  2 (6Y2-2X2-C2Y2R2)] 

N2  R3  1 n 2C  R2 

n 

2 

- -Mf  H(1)(CR)  [X  + — 4—  (3Y2  - X2)]}  (A. 6) 

R ° 2 C R 
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An  error  bound  to  the  first  order  approximation  to  P(a;X;Y)  can 


be  obtained  as  follows: 


00 

~^~2  [ (-u.-iC  )exp(-u  X-iXY)dA| 

tin  ioo  *■  n A 


- ~ V 

tt  In2  I 


r k 

J ]u  _i£  | exp[-X  Refu^JdX 


[Re(C  ) < 0], 


((Re  C2)*  [ReU2>  > 0] , 


then  it  can  be  shown  that 


Re^)  - Re(A2-C2)i4  > j 


(X  < <5), 


VjX^fi2)1*  (X>6). 

Now  <X2-62)**  > X - 5 , if  X i 5.  So  it  follows  that 


Re(Ul)  > 


0 (X  < 6), 
X-<5  (X  > 6). 


Now  U2  and  -i£n  are  in  the  same  quadrant  for  all  real  values  of  X , 
so  j u2~ i Cn | ^ It  | • Using  these  relations,  (A. 7)  becomes 


|£pJ  i 


[ X3dX+  / X3exp l-X(X-d) ]dX 


n N ? ’ J K 

1 n 1 v o 6 

4[2+2$X  * iV  ± 63X3/3) 

N2cJx3 


Similarly  the  error  in  the  second  order  approximation  to  P(a;X;Y)  is 
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ISil  ' I 


iTTN  o 


» 2 

f (-u  -U,  + )exp(-u  X-iXY)dA  | 

i.  n zi,  i. 


r 


Tt|n  I o |2C  (u  -ic  ) 

n i n 


2_  exp  (-XRc  ) dA 


{ J X^dX  + J X^exp [— X ( X— 6 ) ] dX} 


tt|n2C3|  o 

1 n 

2 [24  + 2 46X  + 1262X2  + 463X3  + gV  + 65X5/5] 

IT  I N2C3  | X5 
n 

By  comparing  (A. 10)  with  (A. 11)  , one  can  see  that  the  error  in  P 

2 2 

is  decreased  by  adding  the  second  term  if  | n | X > 6 . 


(A. 11) 
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APPENDIX  B 


In  this  Appendix  an  expression  for  Wo(a;Y)  is  derived.  WQ(a;Y)  as 
given  in  (26)  can  be  manipulated  as  follows 


w f_.Y1  = f e*P.(-iXY,)dX__ 

ota’  ' J UjtUj-i/n) 


/exp(-iAY)dA 

.2  ,2 


♦ i/ft  fex PC;iAV-2dA_ 

-y 


-oo  X - Ap 


W0i(a;Y)  ♦ i/n  W02(a;Y) 


(B.  1) 


•^(aiY)  is  found  by  deforming  the  contour  into  the  lower  half  plane  and 
evaluating  the  residue  to  obtain 

Va;Y)  ' 


-y1  exp(iA  y) 
P P 


(B.  2) 


where  O^arg  Ap<ir  . On  the  other  hand  we  can  rearrange  WQ2  to  obtain 

w fexp(-iAY)dX  _ f°  exp.(-iAyJdX_  , 

W02(a>Y)  2An  W u,  (A-A  ) J u,  (A+A  ) f 
P -00  1 P -00  1 P 


{ijCalY)  - I2C<x;y  ) } 
p 


The  integral  I.  can  be  written  as 


(B.3) 


r^expt-iY(X-Ap)]-l)  / dx  , 

Ij(a;Y)  - exp(-iApY){J^  Ul(A-Xp)  d*  Xp7  } 


««p(-iy )( -i|wc**v  £ * xp  1 
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Y 

IjCa.Y)  = exp(-iA  y){tt  f exp(isA  )hJ1}  fcs)ds  ♦ A f 2 } (B.4) 


y-co  u,(A2-A2) 
1 P 


Similarly  can  be  written  as 


I2(a,Y)  = exp(iApY)ju  / exp(-isAp)H^  (£s)ds  - Ap  J 2 } (B.5) 

-oo  U1  ^ ”^p^ 

By  replacing  these  expressions  for  1^  and  I2  into  (B.3),  we  obtain 

Y (D, 


WQ2(a;Y)  = (ViAp)  J sin[Ap(Y-s)  ]hJ1J  (Cs)ds 


♦ cos(ApY)W02(a;0) 


(B.6) 


where 


"02<a;0)  * / 


dA 


(A2-A2) 

1 p 


An  expression  for  WQ2(a,0)  was  found  by  Olsen  and  Chang  [20]  in 
which  the  contour  is  deformed  around  the  branch  cut,  and  to  this  integration 
was  added  the  residue  of  the  pole  at  Ap.  The  resulting  expression  is 

WQ2(o;0)  * <ft/iAp){2[*n(0-£n(l/fi-iAp)]  * 1*}  (B.7) 


The  principal  value  of  the  logarithms  are  chosen.  The  substitution  of 
(B.2),  (B.6),  and  (B.7)  into  (B.l)  then  yields 

W (a;y)  * (2/A  )cos(A  Y)  [ln(c)  - tn(l/n  - iA  )]  - (tt/A  )sin(A  Y) 
o P P P P P 

♦ »»3(a,Y)  (B.8) 

where 

W,(a;Y)  - (w/x  S)  [ sinlAtt-snH^OsJds 
■*  P P 0 


(B.  9) 
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This  can  also  be  written  as 


W (a;Y)  = (TT1/2A  n){exp(-iA  Y)W,U,Y,iA  ) 

J P p A D 


where 


- exp(iApY)W4(i;,Y,-iAp)} 


Y 

r 


W4(C.Y,t)  - j exp(ts)HQ(Cs)ds 


(B 


(B 


This  integral  can  be  evaluated  by  expanding  the  exponential  into  a power 
series  and  integrating  the  series  term  by  term.  The  resulting  expression 
is 


W4(e,Y,t)  = E 


m=0 


(tYA  j 

m!  m 


(a) 


where 


I (£Y)  = i smH(1)(^Ys)ds 

ul  -j  O 


(B 


(B 


These  integrals  can  be  found  from  the  recursion  relation  for  m > 2 


I^UY)  « (CY)  AHjA,(CY)  + (m-l)(CY)~V1)(5Y)  - (m-l)‘ (CY)'^I  (a) 


,-2„ 


m-2 


IQ  and  1^  can  be  expressed  in  closed  form  as 

I <CT)  - 1}  (?Y)  + (tt/2)  [S  (?Y)H.(1)(CY)  - S (£Y)H(1)  (£Y)  ] 

o o o 1 1 o 

Here  (x)  is  the  Struve  function  of  order  j , and 

Ij/OO  = (a)'1[Hj1)(CY)  + 2i/ (7r?Y)  ] 

These  results  were  obtained  from  Olsen  and  Chang  [20]. 


(B 


(B 


(B. 


.10) 

.11) 

.12) 

.13) 

.14) 

.15) 

.16) 
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APPENDIX  C 

In  this  Appendix  an  expression  for  W (a;X;Y)  is  derived.  W (a;X;Y) 

X X 

as  given  in  (26)  can  be  written  as 

X 00 

W^(a;X;Y)  - - j"  exp(is/n)  j [expf-su^-iAY)/^]  dAds 
o '°° 

X 

- -in  J exp(is/n)H^^  [£ (s^+Y^^lds  (C.l) 

* o 

In  the  case  that  Y ■ 0,  this  can  be  expressed  as 
Wx(a;X;0)  - -inW^S ,X, i/n) 

where  an  expansion  for  is  given  in  Appendix  B,  in  equation  (B.12)  through 
(B. 16) . 

An  expansion  for  is  the  region  Y large  compared  to  X is  found  as 
follows.  Consider  the  function  f (0)  defined  by  f(0)  = ' (£0^) . This  can 

be  expanded  into  a Taylor  series  around  0^  as 

- f^e  ) 

f (0)  - Z r2-  (0-0  ) (C.2) 

m-0  m‘  C 

The  derivatives  of  f with  respect  to  0,  are  given  by  the  recursion  formula 
for  m _>  2 

f tm]  (0)  - -0_1[(m-l)f^"^(0)  + (C2/4)£"“*i(0)J 

where 


f ' (0)  - (-c/20ls)Hj1)(cels) 

Now,  the  substitution  of  0 ■ s^  + and  0£  • X^/2  + Y^  yields 
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By  inserting  (C.3)  into  (C.l)  we  have 

00  Uo  2 / ^ + y 2^ 

m! 


w fa;x,Y)  = -iir  l + X2m+1  I (ix/n) 


(C.  4) 


m=0 


where 


1(0 

m 


= jV- 


1/2)  exp(ts)ds 


(C . 5) 


The  integrals  Im(t)  can  be  expressed  in  terms  of  known  functions  through 
the  use  of  the  recursion  formula  for  m^.2 


■1,  -m 


a t {2  [exp(t)  -(-1)  ] - (2m/t)[2  mexp(t)] 


with 


- (2m-l)I  (t)  - (m-l)I  ,(t)]} 

m-1  m-2 


IQ(t)  = [exp (t)  - l]/t 


(C.6) 


Ix (t ) = (exp(t)  - (2/t)[exp(t)  - IQ(t)]}/t  - IQ(t)/2 


However  for  large  m and  small  t this  method  leads  to  a large 
amount  of  roundoff  error  because  of  the  large  number  of  cancellations.  A 
better  method  for  calculating  I for  large  m and  small  t is  to  expand 
the  exponential  into  a Taylor  series  around  zero  to  obtain 

, k 


'.<«>  • l 


k=l 


(t)' 
k!  Vk 


(C.7) 


where 


K . * f Cs2-l/2)mskds 
m,k  Jq 


which  can  be  computed  from  the  recursion  formula  for  m>l 


2-m  - mK 


m-1  ,k 


m,k  k + 2m  + 1 


(C.8) 
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with 

K . = l/(k+l) 
o,k 

We  note  that  I (t)  is  not  a function  of  a or  Y,  so  these  can  be 
m 

computed  once  for  finding  several  roots  with  constant  n and  X. 

We  also  note  that  the  series  for  W^  given  in  (C.4)  converges 
similarly  to  the  series 


m=0  \ X2  + 2Y2  / 

This  series  converges  for  all  Y > 0,  even  though  the  rate  of  convergence 
decreases  for  small  Y. 
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APPENDIX  D 


In  this  Appendix  an  expression  for  W(a;X;Y)  is  found  that  is  good 
in  the  region  Y smaller  than  X.  The  cosine  term  in  (31)  is  expanded 
into  a power  series  of  Y,  and  this  series  is  integrated  term  by  term 
to  obtain 


W(a;X;Y)  = E 
m=0 


(-1A 


(2m)!  Im(Cl;X) 


where 


A2mexp (-u. X) 


//\  expv-u  a; 

— 7 d\ 

u. (u  -i/n) 

- CO  1 1 


where 


[32/3X2+£2]m  W(a;X,0) 


T exp (-u. X) 

“(a;X-0)  ■ J aT(u  -lVn) 

oo  1 1 


An  expression  for  W(a;X,0)  is  found  by  setting  Y = 0 in  (27)  and 

(29),  and  inserting  these  into  (26)  to  obtain 

X 

W(a;X,0)  - exp (— iX/n)  {- iir  j exp(is/n)H^^  (Cs)ds 


+ (2/A  ) [£n(0  - ?.n(l/n-iA  )]}  (D.2) 

P P 

A series  expansion  for  the  finite  integral  above  is  given  in  Appendix 

B,  in  equations  (B.ll)  through  (B.16).  Let  I be  written  as 

m 

I (a;X)  = A2mW(a;X,0)  - IttR  (a;X)  (D.3) 

m p m 

A recursion  relation  for  R is  found  to  be  for  m > 1 

m — 
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' 


R (a;X)  = (A2  + £2)R  .(a;X)  + f . (a;X)  (D.4) 

m p m- 1 tn- 1 

where 

f (a;X)  = [32m/3X2m] [3/3X  - i/h]H(1)(£X) 
m o 

and 

R (a;X)  = 0 
o 

Inserting  (D.3)  into  (D.l)  one  obtains 

00 

W(a;X;Y)  = cos(A  Y)W(a;X,0)  - itr  E 
P m= 

We  note  that  for  small  £X  the  function  f is  asymptotic  to 

m 

(~l)m(2m)  ! / (X2m)H^  (^X) . This  implies  that  the  series  in  (D.5)  converges 
similarly  to  the  series 

oo  9m 

v (Y/xrm 

, (2m) (2m-l) 
m=l 

This  series  converges  if  and  only  if  Y < X,  so  we  can  infer  that  the  series 
in  (D.5)  converges  in  the  region  Y < X. 


(-iA 


2m 


(2m)! 


R (a;X) 
m 


(D.5) 


. 
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APPENDIX  E 

In  this  Appendix  higher  order  terms  in  the  approximation  of  Q will 
be  derived  and  error  bounds  for  the  first  few  terms  will  be  found.  As 
with  P an  (I4+-l)th  order  approximation  to  Q is  found  by  extracting  the 
first  M+l  terms  in  the  expansion  of  U2  in  (24)  from  u^  in  (23).  Q is 
then  written  as 


M 

Q(ot;X,Y)  = Z Qk<a;X,Y)  + eQM 


(E.l) 


where  Q,  is  the  portion  due  to  the  kth  term  in  the  expansion  of  u„.  Q 

K Z C 

is  given  in  (25),  and  for  k > 0 


2 2 

-2ia  n C 


Qk(a;X,Y)  = 


11  V-».  /A  k 

rrrt  (h) 

(n  -l)nk!  'n  ' 


2k 


I (a;X,Y) 
k 


(E.  2) 


where 


and 


s ■ (l><- 


Ik(a;X,Y)  = (iir)-1  j (A2-Ap)k“1exp(-u1X-iAY)dX 


. 00 

(iTT)"1{[32/9X2+n"2]k"1[-8/ax]  / [exp(-u  X-iAY)/u  ]dA> 

/ nr* 


{ 02/8X2+n~2]k_1i;X(X2+Y2)",SH^1)  [C (X2-«-Y2) } 


In  particular,  Q is  found  to  be 
1 * 

2/s 

Q («;X,Y)  = 9-  H(1)(CR) 

1 (n  -l)n  R L 


(E . 3) 
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The  error  in  approximating  Q by  is  given  by 


-in^/n-u. 


QO1 


?n  r 1 ' l,_u2 

-|  — 4 I 2 — 2~~  exp(-u  X-iXY)dX| 

(n  -1)  -<*>  A -A  1 

P 


<2i  ° 


- V ' 4 

(n  -1) 


/. 


00  exp  [-XRe(u^)  ] 
|u2~in2/n| 


dX 


(E.4) 


2 2 2 a 

If  I (a  ) < I (n  ) , then  both  u„  and  -in  /n  are  in  the  fourth  quadrant 
m — m 2 

2 2 

in  the  complex  plane.  This  implies  that  |u2~in  /n|  > |n  /n | . Using  the 

bounds  on  Re(u  ) given  in  Appendix  A,  the  error  bound  on  Q becomes 
1 o 


I < * I a n I (1+&X) 

"*>■  - ’V-i>„2'  x 


(E.5) 


where 


'(Re,2)54 


(&  : > o), 

(R<£2  < 0). 


Similarly  the  error  in  approximating  Q by  Qq+Q.  is  bounded  by 


- 2a 

i i . 2 1 an 

|eq,' 


OO  |X2-X2| 

;|  I j — "j  exp[-XRe(u  )]dX 


C1  " (n  -l)n  lu^-in  /n|4 

Using  the  same  inequalities  used  in  obtaining  (E.5),  then  (E.6) 
reduces  to 

- 2a3  [2+26X+(52X2+<S3X3/3  +| X2 | X2(6x+1) ] 

|E  I < 2|_aj! 1 E 

V " J(n4-l)n6'  X3 


(E.6) 


(E.  7) 


By  comparing  (E.7)  with  (E.5)  it  is  evident  that  the  error  is  decreased 
2 2 

by  including  if  | n | X > 1 . 
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